Black hole collapse simulated by vacuum fluctuations with moving semi-transparent mirror 
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Creation of scalar massless particles in two-dimensional Minkowski space-time — as predicted by the dynami- 
cal Casimir effect — is studied for the case of a semitransparent mirror initially at rest, then accelerating for some 
finite time, along a trajectory that simulates a black hole collapse (defined by Walker, and Carlitz and Willey), 
and finally moving with constant velocity. When the reflection and transmission coefficients are those in the 
model proposed by Barton, Calogeracos, and Nicolaevici [r(w) = —ia/(u> + ia) and s(w) = uj/(lu + id), 
with a > 0], the Bogoliubov coefficients on the back side of the mirror can be computed exactly. This al- 
lows us to prove that, when a is very large (case of an ideal, perfectly reflecting mirror) a thermal emission of 
scalar massless particles obeying Bose-Einstein statistics is radiated from the mirror (a black body radiation), 
in accordance with results previously obtained in the literature. However, when a is finite (semitransparent 
mirror, a physically realistic situation) the striking result is obtained that the thermal emission of scalar massless 
particles obeys Fermi-Dirac statistics. We also show here that the reverse change of statistics takes place in a 
bidimensional fermionic model for massless particles, namely that the Fermi-Dirac statistics for the completely 
reflecting situation will turn into the Bose-Einstein statistics for a partially reflecting, physical mirror. 

PACS numbers: 03.70.+k, 04.62,+v 



I. INTRODUCTION 

The Davies-Fulling model HUH] describes the creation of scalar massless particles by a moving perfect mirror following 
a prescribed trajectory. This phenomenon is also termed as the dynamical Casimir effect. Recently, the authors of the 
present paper introduced a Hamiltonian formulation in order to address some problems associated with the physical 
description of this effect in the time interval while the mirror is moving JH 01 ; in particular, of the regularization procedure, 
which turns out to be decisive for the correct derivation of physically meaningful quantities. A basic difference with 
previous results was that the motion force derived within the new approach contains a reactive term — proportional to the 
mirror's acceleration. This term is of the essence in order to obtain particles with a positive energy all the time while 
the oscillation of the mirror takes place, and which always satisfy the energy conservation law. Those result followed 
essentially from the introduction of physically realistic conditions, e.g. a semi-transparent or partially transmitting mirror, 
which is perfectly reflecting for low frequencies but becomes transparent to very high ones. 

Here we will study a different aspect of the introduction of physically plausible, semitransparent mirrors, namely the 
particle spectrum produced — in the conditions of the Fulling-Davies effect — by a mirror of this sort which is initially 
at rest, then accelerates during a large enough (but finite) time span, uq, along a trajectory that simulates a black hole 
collapse, as defined by Walker [5], and Carlitz and Willey ]6[]: 

^^(l-e- fcl1 ) (1) 

(in light-like coordinates, where k is some frequency), and finally, for u > uq, is left alone moving with constant velocity 
in an inertial trajectory. 
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We will be interested in calculating the radiation emitted by the mirror from its back (e.g. right) side. As is well-known, 
a perfect mirror that follows this kind of trajectory produces a thermal emission of scalar massless particles obeying Bose- 
Einstein statistics. More precisely, for 1 < tv'/k < e kuo and 1 < lj'/u < e ku ° , one has QHIl] 
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Turning to the case of a partially reflecting mirror — in which we will be mainly interested in this paper — in order to obtain 
the radiation on its right hand side (rhs), we also need to calculate the Bogoliubov coefficient: /3 = 4> % J) L )* . 

We thus first obtain the 'in' modes on the rhs of the mirror when the reflection and transmission coefficients are 
"* and s(w) — , with a > 0, that is, when the Lagrangian density is given 
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where z = g(t) is the trajectory in the (t, z) coordinates. 
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A. The main results 



We now state the main results that will be obtained in this paper (for 1 < uj'/k < e ku ° and 1 < uj'/uj < e ku °). 
1. For a perfectly reflecting mirror, i.e., when u>' <C a, we will see that 
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namely, a thermal radiation of massless particles obeying Bose-Einstein statistics is produced. 
2. For a perfectly transparent mirror, i.e., when a = 0, we will get that 



l/& 2 = o. 



(5) 



In other words, there is no particle production. 

3. In the physically more realistic case of a partially transmitting mirror (transparent to high enough frequencies J^), 
i.e., when a <C w', what we obtain is 
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And, since 
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, we then conclude that a semitransparent mirror emits a thermal radiation of scalar 
massless particles obeying Fermi-Dirac statistics. 

Moreover, Eq. (O does show that there are no ultraviolet divergences when the mirror is semi-transparent (in 
agreement with previous conclusions in 0, Q) and, consequently, the number of produced particles in the u> 
mode, namely J\f u , is finite. In fact, with good approximation we have calculated that 



AC 



2-kuj (fc) (' 



2-kuj /k 



1 



(7) 



It must be here remarked that this phenomenon will in no way occur for a perfectly reflecting mirror, where the 
number of produced particles in a prescribed mode diverges linearly with time (Hill]. 

Actually, the same kind of effect — but now reversed — occurs if a bidimensional/erm/on/c model of massless particles 
is considered, where in the perfectly reflecting case the mirror emits a thermal radiation of fermions obeying the Fermi- 
Dirac statistics. When the mirror becomes semi-transparent (the physically realistic case), the emitted thermal radiation 
will obey Bose-Einstein statistics. We interpret these results as a proof of the fact that the spectrum of the radiation 
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produced by a mirror which follows a trajectory that simulates black hole collapse does not depend on the statistics of the 
field, being just determined by the interaction of the mirror with the radiation field. In our case this interaction is given by 
the reflection and transmission coefficients which depend on the parameter a, that determines the spectrum of the emitted 
radiation. 

Here it is important to emphasize that the word 'statistics' refers all the time to the /3-Bogoliubov coefficient charac- 
terizing the spectrum of the radiated particles and not to the algebra obeyed by the creation and annihilation operators, 
that always satisfy the corresponding canonical commutation relations (anti-commutation, in the fermionic case). That is, 
e.g. in the second case studied the original particles continue to be fermions, but the spectrum of the radiated emission 
corresponds to bosonic ones, when the mirror becomes physical, that is, semi-transparent. As a consequence, fermionic 
number is not actually violated. 

Related phenomenons of similar kind have been reported to occur also in other situations: 

1. In the case of an electric charge following the trajectory v = j(l — e~ ku ). When the radiation field has spin 1, the 
radiation emitted by the charge obeys Bose-Einstein statistics, but when a scalar charge, and consequently an scalar 
radiation field, is considered, the emitted radiation will obey Fermi-Dirac statistics (see d 1 411 for more details). 

2. When measuring the spectrum of a scalar field by using a DeWitt detector which follows a uniformly accelerated 
world-line in Minkowski space-time, one can show that, when the dimension of the space-time is even the Bose- 
Einstein statistics is obtained; however when this dimension is odd the reverse change of statistics occurs (see lfl5ll 
for further details). 

Finally, in an Appendix we will specify some sufficient conditions in order to ensure the convergence of the total 
number of produced particles and of their associated energy. 

II. PERFECTLY REFLECTING, MOVING MIRROR 

Consider a massless scalar field <fi in two-dimensional Minkowski space-time interacting with a moving mirror. Assume 
that the mirror trajectory simulates a black hole collapse HKkI, this is, that it reduces to the following form in the light- 
like coordinates u = t — z and v = t + z: 

{u if u < 

i(l - e~ ku ) if < u < u (8) 
V(u ) + A(u - u ) if u>u , 
with A — e~ kua , where k is a frequency and uq ^> 1. Note that this trajectory can also be written as follows 

{v if v < 

-|ln(l - kv) if < v < v (9) 
U(v ) +A- 1 (v~ v ) if v > v a . 
For a perfectly reflecting mirror, the sets of 'in' and 'out' mode functions are lfl7ll 



f w) = ^3= (e~ luJV ~ e"^ y (")) 9(v - V{u)) 



(10) 



and 



j>™* R [u,v) = -^=L= ( e -^ u - e - iuU ^) 6{v ~ V{u)) 

(11) 

, ^ (u ' v) = -^tm " e~^ y(u) ) 0(u - U(v)), 
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respectively. Our main objective in this section will be to calculate the beta Bogoliubov coefficient 



C? = (0 with w,o/>o, (12) 

where the brackets on the rhs denote the usual product for scalar fields (see, e.g., Birrell-Davies ifl^l ). In order to compute 
this coefficient we choose the right null future infinity domain, J^. There, we have 
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Taking into account that 1 <C w' /k <C A~ x and 1 <C ui' /k <C A^ 1 , we get with good approximation 
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To obtain an explicit expression for the second term on the rhs, we consider the domain 

D = {z e C|Rez G [0, 1 - A],lmz G [-£,0], with k/w' < e < 1} 
and going through the same steps as in |9J], we easily obtain 
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As a consequence, using that |T (1 + iuj/k) | 2 = sin h('Ju/k) (see 111811 ) we get the announced result, that for a perfect 
reflecting mirror the relevant Bogoliubov coefficient is 
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III. PARTIALLY REFLECTING, MOVING MIRROR 

Now we start by reducing the problem to co-moving coordinates (r, p), that is, those for which the mirror remains at 
rest, r being the proper time of the mirror, and we take p such that its trajectory be given by p = 0. Introducing the 
light-like coordinates (u,v), defined as 

u = r — p; v = r + p, (18) 

we will calculate the mirror's trajectory in the coordinates (u,v). Along this trajectory, the length element obeys the 
identity d] 

dr 2 = du 2 = dv 2 = V'(u)du 2 = U'{v)dv 2 . (19) 
An easy calculation yields then the relations 

u if u < 



v — u(u) — < 



and 



u — v(v) 



f(l - e - ku / 2 ) if < u < u 
u(uo) + VA(u — Mo) if u > uq, 

v if v < 



(20) 



f (1 - Vl - kv) if0<u<u 
v(vo) + A~ x l 2 {v — uq) if v > vq. 



(21) 
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When the semitransparent mirror is at rest, scattering is described by the analytical S'-matrix (see lH, l20ll for full details) 



s(u>) r(u>)e 
r(iS)e 2lujL s(uj) 



-liujh 



(22) 



where x = L is the position of the mirror. The S matrix is taken to be real in the temporal domain, causal, unitary, and the 
identity at high frequencies [20], being r(u>) and s(ui) the reflection and transmission functions, which are analytic and 
such that tend in modulus to —1 and 0, respectively, as u> — » (to and 1, as to — > oo). The 'in' modes in the coordinates 
(u, v) are ifToll 
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9u,R 
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Note that the 'in' modes in the coordinates (u, v), namely <f> m , are defined in the right null past infinity by 



-iujv rin 



and, in the left null past infinity J L , by 



= 0, 
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From these definitions it is clear that g™ k {u, v) = g™f.(u(u), v(v)) with k = R,L are not such modes. However, the 
modes g l " k do actually constitute an orthonormal basis of the space of solutions of our problem. As a consequence, if we 
use the fact that g 1 ™^ k — g™* k , we can obtain the following relation 



(27) 



C; fe = J <kj' X (w')(g% <k ; <t>™ k )g$, k , k = R,L 
with x(cj') the sign function. To be remarked is that Eq. < f27b is to be interpreted as follows 

C^= lim f dw' X (u')(ffi, k l<P™ k )g%, k F x (uj% 



being F\(u>') a frequency cut-off, as for instance A2+ A ^^ 2 . 

To calculate the 'in' modes explicitly, we have chosen the coefficients r(w) = ~2f a 
this case, on the rhs of the mirror we obtain 
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Note that (as already advanced) in the case of perfect reflection, that is when a — > oo, we get 

V 4?r M 

and when the mirror is transparent, i.e., when a — > 0, we have 



u> uq 



(30) 



(31) 
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We are interested in the particle production on the rhs of the mirror, for this reason we must now obtain, for u,u>' > 
& = «r;Wm)\ and 0*£, = (fi'Ci)* (32) 

In order to calculate these products we better choose the right null infinity j£ , because here the 'out' modes acquire a 
simple form, namely 



We start by calculating (3^, = 2% J R du <\%% d u <j) r */ l R , with the result 
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Now, provided that to' <C a, Eq. ( f34b turns into Eq. ( TBI ). As a consequence, we precisely obtain the same behavior as for 
a perfectly reflecting mirror. However, in the case a <C w', 
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Finally, a simple but rather cumbersome calculation yields the results 



r (1/2 + iw/A) , 
(see iH), we conclude that 
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Note that in the case a <C d we indeed obtain the nice feature that the number of created particles in the u> mode, together 
with the radiated energies, are both^zn;fe quantities when uq — > oo, in perfect agreement with the conclusions in [8]. More 
precisely, for a partially transmitting mirror the number of produced particles in the uo mode 



did 
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is approximately J °° duj' 
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. In order to calculate this quantity, we split the domain [0, oo) into two disjoints sets, 



[0, k) and [k, oo), respectively. In the second domain we can carry out the approximation d36t to obtain 
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In the other domain, assuming that k <C 1, we have d <C 1 and thus for incident waves of very low frequency the mirror 
behaves like a perfect reflector. For this reason we can use the formula ( TT4T > and a simple calculation yields 



dd 
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Thus, since k <C 1, we conclude that the number of produced particles in the u) mode is approximately 
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and the radiated energy £ = / °° dujftwNu is, also with good approximation, 
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This completes the proof of all of the statements above. 

It is appropriate to remark that there is a crucial difference with the case d <C a, where the number of radiated particles 
in the u> mode diverges logarithmically with u — > oo. In this situation the physically relevant quantity is the number of 
created particles in the lj mode per unit time to = \{u$ + V{uq)) = \u$. This dimensionless quantity is finite and its 
value is given by H 



lim -J-AC = -(t 
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Finally, it is also interesting to calculate the detector response function iflal . namely T{d), for an inertial DeWitt 
detector following the trajectory z = 0. This function is given by 



■F(co)= dt dt'e- luj{f - t ' ) {0,in\<j)(0,t)<j)(0,t')\0,in) 



(44) 



where |0, in) denotes the 'in' vacuum state. This function is related to the average number of produced particles in the 
w-mode, through the relation ^ 



T{uj) = -AC 



(45) 



Then, for a partially transmitting mirror, from Eq. d4Tl ). we conclude that 



(46) 



However, for a perfect reflecting mirror, this quantity diverges. In this case the relevant function is the detector response 
function per unit time lfl5ll . namely V(lo) = lim to ^oo j-^{d). Using ((43), one then obtains the Planckian spectrum 



V{d) = i (e 2 ™l k - l) 



(47) 
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IV. CASE OF THE DIRAC FIELD 



In this section we consider the Dirac equation, in 1 + 1 dimensions, for a massless field 
where here the Dirac matrices are (see i2llo 



In the variables (u, v), Eq. ( |48T > is 



1 




1 
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and we conclude that its general solution can be written as 
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Consider now again the trajectory defined in Sect. II. The vector normal to the trajectory is n M = — y==(y'(u), — 1) 

and the current vector is given by j 7 ' = 2(|G| 2 , |F| 2 ). Then, for a perfectly reflecting mirror we must impose that the 
normal component of the current vanishes on the mirror ( 02211 . 12311 '). that is, j^n^ = 0, and thus the condition follows that 



V'(u)\Gf 
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From here we can calculate the corresponding family of 'in and 'out' modes 
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In this case the beta Bogoliubov coefficient is given by 

& = = / WCA)V£ , R , 
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where (VCr)* denotes the transposed of the vector ip° u R - 

Performing the calculation in the right null future infinity domain, J^, we obtain 
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Then, as in [9], we easily get 
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and thus, since |T (1/2 + iuj/k) \ 2 = cosh( -^/ fc ) (see lfl8IO . we finally arrive to the anticipated result that 
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A. Partially reflecting, moving mirror 

We start with the two orthonormal basis 



2(u) 



(59) 



where the functions S and u are defined in Sect. III. The normal component of the mirror's current given by each of these 
two functions are J= and — jr, respectively. 

It is interesting to observe the analogy between the fermionic and the scalar cases. Note that the current for a scalar 
field is 

f = i{(t>d t cj>* - tffdth -<t>d x cj>* + tpdzfi, (60) 
and if we choose the orthonormal basis: 

fa, u (v) = -±=e-™M, = -^ e --«(«), (61) 

V VKUJ V VKUJ 

we obtain the same normal component of the current as for the respective fermionic function. We can establish the 
following analogy 



^3 



(62) 



and using this analogy we can construct the g ln modes in the fermionic case as follows: simply replace in formulas ( [231 
and d24l i the functions <pj tU with the functions ipj u, to obtain 
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As in the scalar case, the modes g m are not the 'in' modes in the coordinates (u, v). However using the fact that the modes 
g m constitute an orthonormal basis of the space of solutions, we can obtain the following expression of the 'in' modes: 

(65) 



(66) 



Jr 

If we choose the coefficients r(w) = and s(w) = f. , with a > 0, on the right side of mirror we have 



where 
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and 



i I i 



u < o 
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Note that in the case of perfect reflection, that is, when a — > oo, we have 



2tt V°, 

and when the mirror is transparent, i.e., when a — > 0, it turns out that 



(67) 



(68) 



To calculate the production of particles on the rhs of the mirror we must obtain (3^' R , = (VCfl*; R )* and 
Pu',u> = We start b y calculating /3^. Choosing the right null infinity region J^, we have 

P*'* = S R (4'™R) t 4'ls*R du > then ' for 1 < w '/ fc < ete ° and 1 < W 7 W < efetl °' we 8 et 



2tt uj' + ia uj + to' 

^o—iu'/k /•! /'(I - V^) / / 
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Integrating by parts with respect to the s variable, we obtain 



b r.,r ^ Q / 



r(l + iw/fc), 



and thus, in the case of a partially reflecting mirror, it turns out that 



1 /a\ 2 
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Finally, a simple but rather cumbersome calculation yields the result we were looking for 
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V. CONCLUSIONS 



In this paper we have studied in detail the creation of scalar massless particles in a two-dimensional Minkowski space- 
time (the Davies-Fulling theory) and specifically for the case of a semi-transparent mirror, which reflects low frequency 
modes but is transparent to high enough frequencies, being the reflection and transmission coefficients analytic functions 
of the frequency. The considered mirror is initially at rest, then accelerates, during some finite time, along a trajectory 
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that simulates a black hole collapse (as defined by Walker [5], and Carlitz and Willey [6]), and finally rests moving with 
constant velocity. 

When the reflection and transmission coefficients are those in the model proposed by Barton, Calogeracos, and Nico- 
laevici ITol[Tll[T2ll . namely r(w) — —ia/(u + ia) and s(w) = ui/(u + ia), with a > 0, the Bogoliubov coefficients on 
the back side of the mirror could be computed exactly. This has allowed us to rigorously prove that, when a is very large 
(the case of an ideal, perfectly reflecting mirror) a thermal emission of scalar massless particles obeying Bose-Einstein 
statistics is radiated from the mirror (a black body radiation), in accordance with previous results in the literature. More- 
over, we have also seen that when a is finite (the case of a semi-transparent mirror, that is, a physically realistic situation) 
the surprising result is obtained that the thermal emission of scalar massless particles obeys Fermi-Dirac statistics. We 
have also shown in detail that the reverse change of statistics takes place in a bidimensional fermionic model for massless 
particles, namely, that the Fermi-Dirac statistics for the completely reflecting situation gives rise to the Bose-Einstein 
statistics for the case of a semi-transparent, physical mirror. 

The results we have obtained are absolutely solid — they do not hang on a perturbative expansion or approximation of 
any sort. The physical reason for this surprising change of statistics may be found in the fact that the form of the spectrum 
is actually determined not through the statistics of the field but rather by the specific trajectory of the mirror and by its 
interaction with the radiation field. The same kind of phenomenon occurs in the case of an electric charge following the 
trajectory v = 1(1 — e~ ku ). When the radiation field has spin 1, the radiation emitted by the charge obeys Bose-Einstein 
statistics, but when a scalar charge, and consequently an scalar radiation field, is considered, the emitted radiation will 
obey Fermi-Dirac statistics [ 14J. 

Another situation where this kind of features occurs, is when one measures the spectrum of an scalar field using a 
DeWitt detector lfl6i [l7fl which follows a uniformly accelerated world-line in Minkowski space-time. In this case when 
the dimension of the space-time is even the Bose-Einstein statistics is obtained. However, when the dimension is odd, 
precisely the reverse change of statistics occurs in the emitted radiation lfl5l l24l l25l l26ll . 
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VI. APPENDIX 



We will here derive in detail the total number of produced particles and their energy. First, we start with the case of a 
perfectly reflecting mirror. Assuming that the mirror's velocity converges fast enough to some constant, when |tt| — > oo, 
we have 



oR,R 

' J uj.uj' 



1 jUJ 

2^V w 7 



duer luJu e- luJ ' v(u) 



(73) 



and, integrating by parts, we get 



oR-R 



2iri 



i lulu' I du 



V"(u) 



^ — iuu^ — iuj V[u) 



(74) 



For simplicity, we will assume that the mirror's acceleration is discontinuous at the point u = a. After another integration 
by parts, we obtain 



ryR.R 



1 /—! 1 

V LULU 

2tt (w + o/V"(a)) 3 



2tt 



du 



*»« e -*"'v(<0(y"(a-) _ V"(a+)) 
V'"{u) 3uj'(V"(u)) 2 



(lu + oj'V'(u)) 3 (uj + oj'V'(u)) 4 



^ — iuju^ — iu'V{u) 



(75) 
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From this expression, if we further assume that the mirror's trajectory is asymptotically inertial, that is V' (u) > 0, Vu £ 



(see for example (^), it also follows that 



R,R 



and tuo 



R,R 



are integrable functions in the domain [0, oo) 2 \ [0, l] 2 



Now, we are interested in the production of particles in the infrared domain, that is, we want to calculate 
[0, l] 2 . We write the Bogoliubov coefficient as follows 

R . R \_ nj_ 

with B > 0. After integration by parts, we obtain 

1 V W 



R.R 



^ ue -i(^+B^')u e ~iu'(V(u)-Bu) 



ryR^R 



2ir u) + Buj' 



du{V'{u) - B)e 



o'V(u) 



(76) 



(77) 



and thus, if the function | V'(u) — B\ is integrable in R for some B > 0, it can be deduced that 
are integrable functions in the domain [0, l] 2 . 



A 



R,R 



and huj 



oR,R 



An example of this kind of trajectories is 

V(u) = 



Bu u<0 
V(u) 0<u<u o 
V(uq) + B(u — uq) u > uq. 



However, if we are only interested in the convergence of the function Hui 
trajectories that satisfy 



R.R 



(78) 



in the domain [0, l] 2 , we only need 



f du\V'(u) - Si | < oo and f 

J-oc JO 



du\V'(u) -B 2 \<oo 



(79) 



for some non-negative constants Bi and B 2 (here it is important to remark that one of these constants can be zero, that is, 
it is not worthwhile that the trajectory be asymptotically inertial). To prove that statement, we write 



and assume, for simplicity, that V(0) = 0. After integration by parts, we get the expression 

RR 1 , ; Bi- B 2 1 Vww 



(80) 



(I 



2tt 
1 



'0JU1 



2TTUJ + BoUJ 1 



(uj + _Bi u}')(oj + B 2 ui') 2irui + B x w' 



du{V'{u) - B 2 )e- luu e- w ' vi - u \ 



du{V\u) - B 1 )e _iwu e- iu ' v(u) 



(81) 



which already proves the assertion. 

In conclusion, we have here demonstrated that, for asymptotically inertial trajectories with continuous velocity, the 
radiated energy is indeed finite. However it is also possible that an infinite production of particles with very low frequency 
could take place (an infrared divergence). To remove this divergence we must just assume that the initial and the final 
mirror velocities are the same. 

For completeness, we should comment on the very interesting process of particle creation, for the case of a partially 
transmitting mirrors. In this situation, at high enough frequency the mirror behaves as transparent, and then there is no 
particle production, with independence of the mirror's trajectory. On the other hand, at very low frequencies the mirror 
behaves as a perfect reflector, and then we have the same kind of infrared problems as for the perfectly reflecting case. 
As a consequence, if we are only interested in the case when the radiated energy is finite, we must restrict ourselves 
to consider trajectories with a continuous velocity V'(u), Vu G R, which fulfill the condition d79l ) as, for instance, the 
non-asymptotically inertial trajectory: 



if u < 



V(u) 



id 



- ku 



(82) 



:u ) if u > 0. 
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This case can be discussed along the same lines above. 
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